Dynamics of short time— scale energy relaxation of optical excitations due to 
electron— electron scattering in the presence of arbitrary disorder 
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A non-equilibrium occupation distribution relaxes towards the Fermi-Dirac distribution due to 
electron-electron scattering even in finite Fermi systems. The dynamic evolution of this thermal- 
ization process assumed to result from an optical excitation is investigated numerically by solving 
a Boltzmann equation for the carrier populations using a one-dimensional disordered system. We 
focus on the short time-scale behavior. The logarithmically long time-scale associated with the 
glassy behavior of interacting electrons in disordered systems is not treated in our investigation. For 
weak disorder and short range interaction we recover the expected result that disorder enhances the 
relaxation rate as compared to the case without disorder. For sufficiently strong disorder, however, 
we find an opposite trend due to the reduction of scattering probabilities originating from the strong 
localization of the single-particle states. Long-range interaction in this regime produces a similar 
effect. The relaxation rate is found to scale with the interaction strength, however, the interplay 
between the implicit and the explicit character of the interaction produces an anomalous exponent. 

PACS numbers: 71.23.-k, 72.15.Rn, 71.10.-w 72.15.Lh 



The interplay of strong disorder and electron-electron 
interaction is one of the major issues of contemporary 
condensed matter physics. That problem may well be 
behind the insulator metal transition in two dimensions 
or behind the unexpectedly large persistent current 
observed in experiments as compared to theoretical pre- 
dictions 0- Similarly this interplay is responsible for 
the glassy behavior of the electrons recently investigated 
both experimentally |3| and theoretically |4| . Since theo- 
retically the nonperturbative treatment of both disorder 
and interaction is still a very demanding task, numerical 
simulations may yield important insight into the prob- 
lem. In the present paper we present results of such a 
numerical simulation. 

For photoexcited ordered semiconductors it is known 
that Coulomb scattering is a rapid process |5|. In the 
presence of weak disorder, i.e in dirty metals in the diffu- 
sive regime such process may become even faster |(j be- 
cause the particles diffusively can spend more time close 
to each other, which results in an enhanced probability 
of scattering. This enhancement can in other words be 
attributed to the absence of fc-vector selection rules in 
the scattering process. 

Very little is known about the Coulomb scattering for 
the case of strong disorder. In this paper we show that 
with increasing disorder the localization length of the 
single-particle states reduces drastically and hence the 
scattering probabilities as well. Our result is obtained 
from numerical investigations of the dynamical energy re- 
laxation due to electron-electron scattering in a system 
modeling a disordered metal in the localized regime. Our 
model is related to the quantum Coulomb-glass model in- 
troduced and studied in detail in Refs. Q. In those and 
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FIG. 1: Typical time evolution of the root-mean-squared 
deviation of an initially non-equilibrium occupation from 
the Fermi-Dirac distribution. For the case of long-range 
Coulomb interaction. The inset shows the n a distribution 
at different time instants. Time is measured in units of 
Heisenberg-time, A -1 , energy in units of nearest neighbor 
hopping, J. 



subsequent studies that model has been used to deter- 
mine stationary and equilibrium properties of interact- 
ing electrons in a disordered environment. The interplay 
of strong interaction with strong disorder is also respon- 
sible for the emergence of glassy behavior 0, which 
results in logarithmically slow relaxation processes with 
time scales of the order of a day. 

Here we are interested in the short time-scale relax- 
ation of an initially non-equilibrium occupation number 
distribution which is assumed to be a result of an optical 
excitation. Ultimately we are interested in optical phase 



relaxation due to Coulomb interactions in a strongly dis- 
ordered system. These particular relaxation processes 
will take place around a local minimum of the free en- 
ergy of the Coulomb-glass. They are completed long 
before the system moves from one minimum to a lower 
one. Therefore the logarithmic relaxation times arising 
due to the slow process of finding the global minimum of 
the free energy is out of the scope of the present study. 
The treatment of optical relaxation due to interactions 
is not a trivial subject even in ordered semiconductors 
[f| . In order to gain some insight into these processes we 
here study the first the population relaxation in this time 
regime, typical for processes within a given minimum. 
We are aware, however, that phase and population relax- 
ations are not identical. Nevertheless, the dependence of 
the relaxation rates on disorder and interaction strength 
in the situation envisaged is interesting in itself as far as 
optical phenomena are considered. 

In order to investigate the diffusive and the localized 
regimes as well, we use both Hubbard-type short-range 
interactions and Coulomb-type long-range interactions. 
We already anticipate that the former is more appro- 
priate in the diffusive regime as it roughly incorporates 
the screening effect of the other electrons, although, both 
types of interaction yield qualitatively similar results in 
the localized regime. 

In order to investigate the population relaxation due to 
the particle-particle scattering we consider a simplified 
model of a strongly disordered system described by the 
Hamiltonian that consists of two parts, H = H\ + H2, 
where the single-particle part, H\ reads as 

Hi=J2 e *4 c * + E J ij4 c j C 1 ) 

i i,j 

where Ci (cj) annihilates (creates) an electron on site i, 
(i.e. a state i)). We consider electrons without spin. 
The atomic energy levels, £j, are taken randomly from 
a box distribution of width W around zero mean value. 
Jij describe the hopping amplitude from site i to site j. 
Nearest neighbor approximation has been used, with a 
constant hopping rate J taken as the unit of energy. The 
sites are assembled in a regular one-dimensional lattice 
of unit lattice spacing with periodic boundary conditions. 

The second part of the total Hamiltonian contains 
the two-particle interaction which in site representation 
reads as 

H * = \ E V V ^ ~ K ) ( c fa - K ) ' ( 2 ) 

ij 

where for the sake of charge neutrality we have already 
included a compensating charge of Ke at each lattice site 
where K is the filling factor. The interaction matrix ele- 
ment is either of short-range or long-range type. In the 
former case Vij — Uq when two electrons are on the neigh- 
boring sites, \i — j\ = 1, and zero otherwise. For long- 
range interaction, Vij = Uo/\i — j\, Uq > characterizes 
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FIG. 2: Dimensionless relaxation rate in units of the mean 
level spacing for the case, 7 = F/A, of (a) long-range 
Coulomb interaction and (b) short-range Hubbard interac- 
tion as a function of the dimensionless disorder strength, 
w = W/J. The initial occupation is similar as in the in- 
set of Fig. □ The different curves are labeled according to 
u — Uo/J- The insets show the curves rescaled with the 
dimensionless interaction parameter, u. The system size is 
TV = 20. 

the strength of the repulsion between electrons located at 
neighboring sites. In any case due to the Pauli-principlc 
the electrons are not allowed to occupy the same site. 

The electron-electron scattering is evaluated in an ef- 
fective single-particle basis. This basis is obtained from 
the diagonalization of the Hamiltonian H± including the 
'diagonal' part of the interaction, i.e. as a first step we 
selfconsistently obtain the Hartree-Fock (HF) solution of 
the Hamiltonian Q by replacing the parameters £j and 
Jij as Ei + \ Vij n jj an d J&(i,j) ~ \Vij n ji- I n f ne HF 
basis our original Hamiltonian can be expressed as 

H = ^2e a clc a + u ll c i c l c t c s ( 3 ) 
where the single-particle part, is obviously diagonal in 
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the HF basis with the e a 's being the HF eigenvalues. In 
the residual interaction 

U~a(l — ^ Vij Cja, C*0 Cj 7 CjS , (4) 
ij 

where the {Ci a } numbers represent the HF states in site 
representation. This model is in fact the standard two- 
body model of interacting (TBRIM) fermions. Within 
the TBRIM, for example, the U^t values are chosen ran- 
domly from a Gaussian distribution assuming that the 
single-particle states are sufficiently chaotic, i.e. delo- 
calized. The typical matrix element U ~ A/g, where g 
is the dimensionless conductance of the system and A is 
the mean HF-level spacing. In our case these interac- 
tion matrix elements contain the information about the 
microscopic details of the original model, e.g. the long- 
range correlations due to the Coulomb-interaction and 
the presence of a disorder potential. Hence, especially in 
the case of strong disorder we are not allowed to use the 
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As we have noted already, in the localized regime we can- 
not apply assumptions of ergodic wave functions in order 
to estimate the typical value of U^i. Also our spectra are 
discrete therefore the (^-function in JSJ is not possible to 
be satisfied exactly. However, we may approximate it 
with a box of finite width of the order of the mean level 
spacing, A. This approximation enables us to call the 
relaxation described by Eq. © in fact an inelastic pro- 
cess because the finite width effectively results in loss of 
energy. 

Integration of the above equation using a standard 
fourth order Runge-Kutta procedure gives the time evo- 
lution of n a (t) . The form of a Fermi-Dirac distribution 
it-fd(E) = 1/(1 + exp[(3(E — fi)]) is fitted at every time 
step. This fit provides a 'chemical potential', fj, and an 
'inverse temperature', /3. However, we are more inter- 
ested in the error of this fit, cr 2 (t) = J2 a ( na ~ n FD{£ a )) 2 ■ 
This quantity characterizes how close the distribution n a 
is to npr>. An example for one single realization over 
N = 20 sites with half filling is presented in Fig. ^ where 
we can clearly see an approximately exponential decrease 
that initially characterizes the relaxation process. From 
these initial exponentials a relaxation rate, T, can be ob- 
tained via a(t) Ri do exp(— Tt). 

The scattering probability between the pairs of single- 
particle states {a, (3} and {7, 5} is {U2p) 2 which is explic- 
itly proportional to Uq. Therefore for small enough Uq 



TBRIM model. Note also that @ depends on the inter- 
action strength Uq explicitly through Vij and implicitly 
through the coefficients, {Ci a }. 

This HF basis corresponds to a zero temperature equi- 
librium distribution of the occupation numbers n a = 
(c a c a ) that equals 1 (0) for e a < Ep (e Q > Ep). 

We assume that an excitation process has somehow 
generated an initial, non -equilibrium n a distribution of 
the form n a = Z^ 1 exp [— (e a ~ E c ) 2 /2w 2 ] , where E c is 
the center of the 'excitation' and w is its energy spread. 
Z is fixed by the condition J2 a n a = N. The center of the 
'excitation' is chosen in the lower half of the energy band, 
typically at its middle or at the bandedge. The width, w, 
is typically chosen to be one fourth of the bandwidth with 
which one could more-or-less avoid the possibility of the 
non-physical situation of n a > 1. During the numerical 
simulation this initial distribution is assumed to relax to- 
wards equilibrium via electron-electron scattering. This 
process is described by the Boltzmann equation |6( 



(5) 
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when the Hartree-Fock states differ very little from the 
non-interacting basis we expect T ~ Uq. The power of 
two should, however, be an approximate value, since for 
strong enough interaction and also for strong enough dis- 
order we expect a different exponent due to the implicit 
character of the interaction encoded in the coefficients 

{Cia}- 

The relaxation rate obtained from the exponentials as 
shown in Fig. ^ has been collected and averaged over 
many realizations. In Fig. |2 we show data obtained 
for several interaction strength, Uq, and disorder, W. 
We can clearly identify qualitatively that for the case of 
short-range interactions, weak disorder produces an in- 
crease of the relaxation rate 0] . As disorder is increased, 
however, the rate decreases. Note that a long-range in- 
teraction induces a much faster relaxation as compared 
to a short-range one. 

As we can see in the insets of Fig.|21the relaxation rate 
grows as a power of Uq that is smaller than two. This is 
due to the strong perturbation the interaction makes on 
the HF states as compared to the non-interacting basis. 
One may also detect a slight difference in the exponents 
between the two types of interaction. 

We may summarize that the error of the occupation 
number distribution decays for short times roughly as an 
exponential (Fig. therefore we may expect that the so- 
lution of Eq. (JSJ is also an exponential whose derivative 
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FIG. 3: Average dimensionless relaxation rate, 7 = (r a )/A 
as a function of the dimensionless disorder strength, w = W/J 
as obtained from the exponential ansatz for short times. The 
different curves are labeled according to u = Uo/J. Open 
symbols stand for Hubbard type interaction and filled symbols 
for Coulomb type interaction. The data for the Coulomb 
interaction has been scaled down by a factor of N, where 
N — 20 is the systems size. 

If we substitute this ansatz into Eq. (J5J 
we can calculate the r Q values from the initial deriva- 
tives. This allows obviously for a much better statistics, 
however, will still strongly depend on the initial n a dis- 
tribution. An average over the individual rates, T a , and 
obviously over many realizations of the disordered po- 
tential is presented in Fig. [3J We can observe that the 
effective Uq dependence here is roughly two owing to the 
simple approximation of the exponential ansatz for short 
times. 

In summary our numerical results show that weak dis- 
order indeed causes the energy relaxation of an initially 
non-equilibrium occupation number distribution towards 
a Fermi-Dirac distribution to become faster for the case 
of short-range interactions. On the other hand with 
an increase of disorder the single-particle localization 
volume decreases hence the quasi-particles have smaller 
chance to effectively scatter and therefore the relaxation 
rate decreases considerably. Note that the long-range 
type interaction is more effective, therefore T is orders of 
magnitude larger as compared to the short-range inter- 
action. However, qualitatively both of them produce a 
similar tendency as a function of large enough disorder, 
W. From the results presented here one could expect 
that coherent phenomena at elevated density of particles 
may be more robust with respect to dephasing processes 
in strongly disordered systems. Thus one could suspect 
that these systems may serve as a testing ground of co- 
herent phenomena even at elevated carrier densities, es- 
pecially in the case when these phenomena do require 
the presence of strong disorder anyway. One such phe- 
nomenon is for instance the current echo [8j. 

The loss of coherence, i.e. dephasing, following inter- 



band photoexcitation has been widely studied in semi- 
conductor systems, where the optically induced coher- 
ence is monitored in the time domain by ultrafast non- 
linear optical techniques 0, ^| . Using many-body the- 
ory it has, however, been shown that the dephasing of 
the interband polarization due to the Coulomb many- 
body interaction cannot be estimated directly from the 
carrier-carrier scattering rates but is a rather compli- 
cated phenomenon. So-called in scattering contributions 
largely compensate the dephasing provided by the usual 
out scattering terms 0, H, 0, 0] . As a result the depen- 
dence of the dephasing rate of the interband polarization 
on the carrier density n is rather weak, see Ref. E3 where 
for two- and three-dimensional systems a dependence of 
cx n 1 / 3 was found. It would be very interesting to in- 
clude coherent contributions in our type of approach and 
to investigate to what extent such results are altered in 
the presence of disorder. 
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